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Abstract 

We present the calculation of the order a 2 s contributions to the coeffi- 
cient functions (singlet and non-singlet) belonging to the longitudinal spin 
fragmentation function g\{x, Q 2 ) measured in polarized electron-positron an- 
nihilation. From this calculation we also obtain the two-loop contribution to 
the timelike spin splitting functions P qq and P gq . These splitting functions 
are in agreement with recent results in the literature which were obtained 
by using a different method. We show that in second order the renormaliza- 
tion constant needed for the HVBM prescription of the 75-matrix is process 
dependent. 



1 



1 Introduction 



The measurement of the fragmentation functions F^(x, Q 2 ) in unpolarized 
electron-positron annihilation provides us with a new possibility to test the 
predictions of perturbative quantum chromodynamics (QCD). Here the reac- 
tion under consideration is given by e + + e~ — > H + "X" where H represents 
the detected hadron and "X" stands for any inclusive hadronic state. Be- 
cause of the similarity between the above process and deep inelastic lepton- 
hadron scattering + H — > e~ + "X" one can consider the fragmentation 
functions as the timelike analogues of the well known deep inelastic struc- 
ture functions measured in the latter process. This similarity becomes clear 
if one looks at the experimental and the theoretical characteristics of both 
functions. As far as experiment is concerned one has a full control of the 
background and a large amount of data so that the systematical and statis- 
tical errors are rather small. The data are collected over a wide range of Q 2 
values which enable us to test the scaling violations exhibited by both func- 
tions (for the fragmentation function see e.g. [ij, @, 0. From the theoretical 
viewpoint structure functions and fragmentation functions have the common 
feature that the Born approximation is determined by the electroweak in- 
teractions only. Therefore it is easier to calculate their QCD corrections at 
higher order in the strong coupling constant a s than for other quantitities. 
It turns out that the most complicated virtual corrections to these reactions, 
which are represented by n- point functions with n > 4, appear in much 
higher order than for other processes like e.g. jet production. Furthermore 
the semi inclusive nature of the above processes leads to a considerable sim- 
plification of the phase space integrals. Another important feature is that one 
does not have to deal with problems like jet definitions and hadronization ef- 
fects in the case of inclusive quantities. Because of the these properties it was 
possible to compute the order a 2 contributions to the coefficient functions 
corresponding to the structure functions @, 0, @ and the fragmentation 
functions JF|, M, ||, [TO]. In spite of the similarities between the semi inclu- 
sive quantities in deep inelastic scattering and eletron-positron annihilation 
there is also one striking difference. This can be attributed to the current 
current correlation function which shows up in the expressions for the struc- 
ture functions as well as in the fragmentation functions. In the case of deep 
inelastic scattering the current current correlation function is dominated by 
the light cone in the Bj0rken limit whereas in the case of electron-positron 
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annihilation this property does not exist. Therefore the light cone expan- 
sion techniques cannot be applied to the fragmentation functions. These 
techniques allow us to relate parton densities to the expectation values of 
hadronic operators between hadronic states. It also enables us to compute 
sum rules of structure functions to even higher orders in perturbation theory 
than is possible for the coefficient functions (see e.g. |TT|). 



Up to now we have only discussed unpolarized scattering. It would be 
interesting to see whether QCD also yields good predictions when the hadron 
H is polarized. This has been investigated theoretically as well as experimen- 
tally during the last ten years for polarized deep inelastic scattering. The 
most effort was devoted to the longitudinal spin structure function gi(x, Q 2 ) 
because of the disagreement between the data |T2j and its first moment given 



by the Ellis- Jaffe sum rule ||13|| . This discrepancy lead to a thorough investi- 
gation from which one can conclude that it is more difficult to make QCD pre- 
dictions for spin dependent quantities than for spin averaged ones. Therefore 
we expect that the same will happen for electron-positron annihilation when 
the detected hadron H in the final state is polarized. This process allows us 
to study the timelike analogue of the above longitudinal spin structure func- 
tion which we will denote by the fragmentation function g^(x, Q 2 ). Contrary 
to polarized deep inelastic scattering there do not exist data for the longitu- 
dinal spin fragmentation function. At this moment there are some data on 
polarized A-baryon production (see e.g. fl4}] ) so that in the future one maybe 
able to measure the fragmentation densities AD^(z, /i 2 ) for a = q,q,g which 
appear in gi(x,Q 2 ). On the theoretical side considerable activity has been 
devoted in the past few years to understand the fragmentation in polarized 
A-baryon production in spacelike as well as in timelike processes. For recent 
work we refer to pf, p|, |7|> PI' @ ( for older work see also PI)- A P art 



from the nonperturbative aspects it is also important to calculate the QCD 
corrections to the spin fragmentation functions. Recent progress has been 
made in the computation of the next-to-leading order (NLO) corrections to 
g±(x,Q 2 ). The order a s contributions to the coefficient functions have been 
calculated in JTB|-||2^| and the two- loop contributions to the timelike spin 



splitting functions are presented for the first time in [pOfl . It is our aim to 



extend the above calculations by including the order a s contributions to the 
coefficient functions which is a first step to obtain the next-to-next-to-leading 
(NNLO) corrections to gf(x,Q 2 ). 

Our paper will be organized as follows. In the next section we introduce 
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our notations and present the framework of the computations of the longi- 
tudinal spin timelike coefficient functions. In section 3 we give an outline of 
the calculation. From the latter we also obtain two out of the four splitting 



functions computed in |2(J • Further we discuss the intricacies of the prescrip- 
tion of the 75-matrix in n dimensions and show the process dependence of 
the renormalization constant needed to restore the Ward-identities when this 
matrix is treated in the scheme of HVBM [^T[ , [^] . The long expressions for 
the order a 2 s corrected coefficient functions are presented in appendix A. 
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2 Fragmentation functions in polarized and 
unpolarized electron-positron annihilation 

Single particle inclusive production in electron-positron annihilation is given 
by the reaction 

e-(Ai, h) + e + (A 2 , l 2 ) -> V(q) -> H(«,p) + "X" , (2.1) 

where "X" represents any inclusive final hadronic state. Further Aj, li (i = 
1, 2) denote the spin and momentum of the incoming leptons and s,p stand 
for the spin and momentum of the hadron H detected in the final state. If we 
denote the momentum of the virtual vector boson V, with V — j or V — Z, 
by q = h + h the Bj0rken scaling variable is defined by 



X =%T> q 2 = Q 2 >0, 0<x<l. (2.2) 



The spin averaged differential cross section can be written as 

^ = jEE^ft;^^ > (2-3) 

where V and V run over 7 and Z and iV stands for the colour factor of the lo- 
cal gauge group SU (N) . Further we have denoted the spin averaged leptonic 
tensor by Cfy )V ,Mi] h) and T>(yy^(Q 2 ) stands for the electroweak vector bo- 
son propagators. If we average over the spin of the detected hadron H one 
obtains the hadronic tensor W™ era9e ' n '( y ' v '\p, q) which can be expressed into 

the fragmentation functions F^ V ' V \x,Q 2 ) (see below). An explicit expres- 
sion for Eq. (|2.3|) can be found in and p3| . In the case the electron is 
polarized we can choose Ai =| and average over the spin A2 of the positron 
in Eq. (|2.1| ). The difference between the cross sections where the detected 
hadron is polarized parallel s =| and antiparallel s =f with respect to the 
electron is given by 



dxdfl dxdfl 

N X 

x x '-(V,V) 



T E E £ftvoU> l ^ W vy/) (Q 2 )W^r' n ' {V ' V '\s,p, q) . (2.4) 



V V 



5 



For a polarized electron the leptonic tensor is given by £jyyn(|, h', h) and 
the spin dependent hadronic tensor Wf^ n ^ v ' v '\s,p] q) can be expressed 
into the spin fragmentation functions gf'^ V ' V \x, Q 2 ) (see below). Explicit 



formulae for Eq. ( |2.4|) can be found in | r?| . In general the hadronic tensor 
can be written as 

i-£ J rfV*(0 | J» | n(s,p),X)(K(s,p),X | jf(0) | 0) 

= W^ ra9e ' H ' (V ' V '\p, q) + W°r> U ' (Vr) (s,P, q) • (2-5) 

The spin average hadronic tensor can be expressed into the fragmentation 
functions Fi(x,Q 2 ) as follows 

W™ ra9e > H ' {V > V,) (p, q) = (~9^+ q -f-)Ff^ V \x, Q 2 ) 

+ P ^F^ v \ X ,Q 2 ) +i e^p-F^ v '\ X ,Q 2 ). (2.6) 



p • g 2p • g 

Similarly one can write the spin dependent hadronic tensor as 

w°r' H ' {vy '\s,p, q ) = 

2£/^A CT - 9l {30, Q ) + itfMuXaqxlSa Pa 

v,n ■ ci \ r> ■ a 



2p ■ q " V p ■ q > 2p ■ q 

( P»S V + Sjpy P^Pu\ gf [Vy '\x,Q 2 ) P^py H,(V,V)/ r n 2^ 

-r\ ~ o-q i r s ■ q- — c/ 4 yx, j 

v 2 p • g ' p ■ q \P ■ q) 

+ (_ 9( , + Ss*)Ll» s f.KV) (!Ci0 . )i (2 . 7) 

V f q z / p ■ q 

where <7?(;r, Q 2 ) denote the spin fragmenation functions. Further we have 
introduced the abreviations 

p-q - s-q 

P»=P„- ~^q„ , s M = s M - —q„ . (2.8) 

The above fragmentation functions can be written as a product of the elec- 
troweak couplings, describing the interaction of the vector boson V with 
the quarks, and the fragmentation functions i^ H , gf so that the superscripts 
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(V, V) can be dropped in the latter quantities. The QCD analysis of the frag- 
mentation functions proceeds in the same way as is done for the structure 
functions in deep inelastic lepton-hadron scattering. Hence it follows that the 
leading contributions to the fragmentation functions Ff are all of type twist 
two. The same holds for the spin fragmentation functions gf (i = 1,4,5). 
However gf and gf also receive leading contributions of type twist three. 
This means that the latter cannot be simply expressed into convolutions of 
coefficient functions with fragmentation densities as can be done for twist 
two contributions. Limiting ourselves to twist two one can simply put s = p 
in Eq. ( |2.7| ) so that the functions gf and gf drop out of the equations and 
we are left with the longitudinal spin quantities gf with (i = 1, 4, 5). Further 
one can derive 

Ff(x,Q 2 ) = -Ff(x,Q 2 ) , gf(x,Q 2 ) = -gf(x } Q 2 ) , i = 4,5,(2.9) 

from charge conjugation invariance of the strong interactions. Hence it fol- 
lows that the above fragmentation functions only receive non-singlet (with 
respect to flavour) contributions. The other quantities, given by Ff, Ff and 
gf, also have singlet parts. The QCD corrections to the spin averaged frag- 
mentation functions Ff (i = 1, 2, 3) have been calculated up to second order 
in the strong coupling constant a s in ||, H], ||TU[ (for the longitudinal 
fragmentation function Ff = Ff — IxFf see also @). Up to order a 2 the 
QCD corrections to gf and gf are the same as for the non-singlet parts of 
Ff and Ff respectively. The same holds for the nonsinglet part of gf which 
is equal to the correction to Ff. These equalities will be explained below. 
What is left is the singlet contribution to gf which will be calculated up to 
order a 2 in this paper. This fragmentation function is the only one which 
shows up in the spin dependent differential cross section ( |2.4| ) if V = 7 so 
that process ( [2.1|) is purely electromagnetic. In this case Eq. ( |2.4|) simplifies 
enormously and one gets 

ax a cos ax a cos 6 Q z 

after integration over the azimuthal angle. Notice that gi(x, Q 2 ) can be also 
measured in the reaction where both the electron and positron are unpolar- 
ized. In this case it appears in the cross section (|2.3|) because of the axial 
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vector coupling of the Z-boson to the incoming leptons. In the purely elec- 
tromagnetic case the longitudinal spin fragmentation function can be written 

as 



+ n f ADr S (- z ,Mi)AC?s(z,^)]. (2.11) 



The spin parton fragmentation densities are defined by 

AD*(z, M 2 ) = D*j:(z, M 2 ) - D%(z, M 2 ) , (2.12) 

with a = q,q, g. The quark fragmentation densities can be distinguished into 
singlet (S) and non-singlet (NS) parts. The former is given by 

n f 

ADf s (z, M 2 ) = £ (ADf(z, M 2 ) + ADf(z, M 2 )) . (2.13) 
fc=i 

Here k = 1, 2, • • • n/ means = u,d, ■ ■ -Uf where stands for the heaviest 
light flavour appearing in the sums in Eqs. ( |2.11|) , (|2.13|) . The spin non- 
singlet quark fragmentation densities are denoted by 

AD%' m (z, M 2 ) = ADf(z, M 2 ) + ADf(z, M 2 ) - — ADf' s (z, M 2 ) . (2.14) 

n f 

The fragmentation densities and the coefficient functions ACi >a (z, Q 2 /M 2 ) 
(a = q, g) depend on the partonic Bj0rken scaling variable z (see below) and 
on the factorization scale M which we have put equal to the renormalization 
scale. 

The coefficient functions are computed from the partonic subprocesses 
which can be generally denoted by the reaction 

V(q) -> "a(s,p)" + ai(si,pi) + a 2 (s 2 ,p 2 ) H h ak(s k ,p k ) , (2.15) 

where a represents the detected parton (a = q, q, g) with momentum p which 
fragments into the hadron H. The remaining partons a k belong to the in- 
clusive state and one has to integrate over all momenta p k and to sum over 
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all spins Sk- The reaction in Eq. (|2.15| ) can be described by a partonic 



structure tensor, indicated by a hat, which has the same form as Eq. (|2.5| ) 
where H(s,p) and X are replaced by a(s,p) and {a^} respectively. In the 
case the detected parton in Eq. fl2.15|) is given by a massless (anti-) quark 
the partonic structure tensor for unpolarized scattering can be written as 

W°r a9e ' q ' {V ' V,) (j>, q) = Itf GW(p, q) • (2-16) 

In the above equation denotes the amplitude squared where the Dirac 
spinors of the detected quark q(p) are removed. The indices /i and v refer 
to the (axial-) vector couplings of the vector bosons V and V respectively. 
Further we have integrated over all momenta and summed over the spins of 
the inclusive state in reaction ( 2.15Q . If the gluon is detected in the final 



state of reaction (|2.15 ) we can write 



w;r a9e ' 9 ' (v ' v '\p, q) = PapG$(p, q) • (2-17) 

Here G^£ is defined in an analogous way as G^ u above and a,/3 refer to the 
Lorentz indices of the polarization vectors of the detected gluon i.e. e a (p), 
e^ip) which are removed from G?£. In unpolarized scattering P a p stands 
for the spin averaged sum over the polarizations. The same can be done for 
polarized scattering. Choosing s = p the partonic spin structure tensor reads 

W^ v ' v '\s,p, q) \ s=p = TrrrfG^ip, q) . (2.18) 



If the gluon is detected in the final state of reaction ( 2.15 ) we can write 



W™ era9e ' 9 ' iVr) (P, q) = AiV3#(p, q) . (2.19) 

In polarized scattering AP a/ 3 stands for the difference between the right- and 
the lefthanded polarizations. Here we will choose (see ||) 

AP Q/3 = — £ q/3(7 apV • (2.20) 
p-q 

From Eqs. (^Wj )- (^I% one can infer the partonic fragmentation func- 
tions Fi t a(z,Q 2 ) (i = 1,2,3) and gi, a (z,Q 2 ) {i = 1,4,5) with a = q,g and 
z = 2p ■ q/Q 2 (see Eq. (|2.15|) ) is the partonic Bj0rken scaling variable. They 
are defined in a similar way as the hadronic fragmentation functions in Eqs. 
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[B[), ( |2.7| ). Notice that the partonic fragmentation functions are collinearly 



divergent due to the higher order QCD radiative corrections. These diver- 
gences have to be removed via mass factorization so that one obtains the 
coefficient functions as e.g. shown in Eq. (|2.11|) . The calculation of the spin 
averaged quantities !F i a has been carried out up to order ot 2 s in 0, ||, ||, [|10| 



Here one has applied the method of n-dimensional regularization to compute 
the many body phase space integrals and the loop integrals appearing in G^ v 
(Eq. (|2.16|) ) and G?£ (Eq. (|2.17]) ). We will use the same method to com- 



pute the partonic spin fragmentation function g\ A which will be presented in 
the next section. In this case one has to find a suitable prescription for the 
75-matrix and the Levi-Civita tensor in Eq. ( [2.20|) . Here we have adopted 



the HVBM prescription given by 't Hooft and Veltman [21] which has been 
worked out in more detail by Breitenlohner and Maison |2j| (see also [p5|[ . 
[BBj). The computation of g^ a is completely analogous to its deep inelastic 
counterpart described in M and the spin averaged partonic fragmentation 
functions in f| so that we do not have to repeat the details in this pa- 
per. Before presenting the results in the next section one can derive several 
relations between the coefficient functions in the polarized and unpolarized 
case. However they are only valid up to order a 2 . The non-singlet coefficient 
functions satisfy the following equalities 

A/^NS _ /7NS A^ NS — ^ NS A^ NS — ^ NS (O 01 ^ 

For the singlet parts we have 

AC 4 p * = , AC^ = 

AC 4 s 9 = , AC 5 s , = 0, (2.22) 



with the definition 



ACl = AC^ + AC™ . (2.23) 



The same relations also hold for the partonic fragmentation functions from 
which the coefficient functions are derived via mass factorization. The re- 
lations in Eq. ( |2.21| ) follow from the anticommutativity of the 75-matrix 
which is manifest for a regularization method in four dimensions. In the case 
of n-dimensional regularization it occurs after the introduction of additional 
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renormalization constants which restore the various Ward-identities violated 
by the prescription used for the 75-matrix and the Levi-Civita tensor. The 
partonic fragmentation function gf^ is derived from Eq. ( |2.1&| ) where the 



indices n,u either stand for 7^,7^ or for 7^75, 7^75. Further if originates 
from so called rainbow graphs (see e.g. fig. la in [fLl]| ) then the 75-matrix 
in Eq. (|2.18| ) can be anticommuted in the string of 7-matrices so that either 



the index fi or the index v represents an axial-vector. In this way Eq. ( j2.1§| ) 



becomes equal to Eq. ( |2.16|) which leads to In the case that either fi or 
v (but not both) represent 7^75 or 7^75 one can anticommute the 75-matrix 
in Eq. ( |2.18| ) till it annihilates the other one in G^ v and one obtains Eq. 
( |2.16|) . In this way one obtains gf^ = F^q an d 95 q = -^lf- The relations 
in Eq. ( |2.22| ) follow from charge conjugation invariance of QCD so that one 
gets 

^3,a = ~^3,a , 9i,a = ~9i,a , 1 = 4, 5 . (2.24) 

Notice that the relations in Eq. (|2.21|) do not hold anymore beyond order a 2 . 
A nice example is given by the graphs in fig. lb of |TT| (light by light scatter- 
ing graphs) appearing in deep inelastic lepton-hadron scattering (spacelike 
process). These graphs contribute in order to J-f^ but not to gf^. In first 
order of a x the coefficient functions AC^ and AC^r were calculated in HIT 



and they satisfy the relations in Eq. (|2.21|) . The order a 2 contributions are 
equal to C^f an d respectively and can be found in B. 
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3 Order a 2 s contributions to the coefficient 
functions AC\ a 



In this section we present the results for the coefficient functions ACi jCl which 
emerge from the computation of the partonic fragmentation functions gi ja 
with a = q,g. The latter are collinearly divergent and the singularities are 
regularized via n-dimensional regularization so that they manifest themselves 
as pole terms of the type (l/e) k , with e = n — 4, in gi <a (z, Q 2 , e). For the 
definition of the 75-matrix we use the HVBM prescription [^2j ( see 
also [23], |2D|]). This prescription destroys the anticommutativity of the 75- 
matrix in Eq. ( |2.18| ) so that gf q 7^ q . This will lead to a violation of the 
relations between the coefficients in Eq. ( [2.21| ) and the nonrenormalization 
properties of the non-singlet axial- vector current. Furthermore as is shown in 



lPq| the Adler-Bardeen theorem |27j will be also violated. These properties 
can be restored by introducing an additonal renormalization constant Z r qq 
(r = NS, S) which has the following unrenormalized form 



'qq 



6(1 



"Arc' 



+ & 2 Sl^Y\ - -JozS(z) + zlf\z) - 1^ Q \z)\ , (3.1) 



up to order a 2 s . Further we have defined the spherical factor as 



exp 



\e{ lE 



In 47r) 



(3.2) 



which is an artefact of n-dimensional regularization. The same holds for 
the scale [i which originates from the dimensionality of the gauge coupling 
constant in n dimensions (g — > g(/i)~ £ / 2 ). The mass parameter M denotes 
the renormalization scale which will be put equal to the factorization scale. 
Further (3 is the lowest-order coefficient in the series expansion of the (3- 
function which, up to order g 5 , is given by 



M = -A) 



167T 2 



01- 



167T 



2\2 



00 
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C A --n f T f , g 2 = 47ra s ,(3.3) 



where Ca = N and Tf = 1/2 are colour factors in SU(N). Notice that in 

+ UfZyL 3 '^ where UfZ^l 3 '^ is coming from the 



second order Zg q ^ 



X NS,(2) 
"99 
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NS 



light by light scattering graphs in fig. 1 of [p5 |. In the non-singlet case Z 
can be derived from 



f™{z,Q\^) 



\Q 2 =M 2 



The first order coefficients are given by 



^qq 

with C F = (N 2 - 1)/2N and 



&W = C F [ - 8(1 - z) Inz - 4(1 - z) ln(l - z) + 10 - 8^ 



(3.4) 



(3.5) 



(3.6) 



The second order coefficient z^ q s '^ is given in appendix A. We will now 
present the general form of the partonic fragmentation functions which follow 
from mass factorization and renormalization. The expressions are written in 
such a way that the non-logarithmic parts of the coefficient functions denoted 
by Ci t k (k = q, g) and the the DGLAP splitting functions are presented in 
the MS-scheme. This scheme holds for coupling constant renormalization as 
well as mass factorization. 

In the Born approximation the partonic reaction in Eq. ( 2.15 ) is given 
by (see fig. 2 in M) 



V -*q + q. 



(3.7) 



where either the quark or the anti-quark is detected in the final state. This 
leads to the contribution 



3 (0) 

9l,q 



l = S(l-z). 



(3.8) 



In order a s we have the following contributions. First we have to compute 
the one-loop virtual corrections to the Born reaction (|3.7p (see fig. 3 in 0) 
and add them to the gluon bremsstrahlung process (see fig. 4 in M) 



V -> q + q + g . 



(3.9) 



In the case the quark or the anti-quark is detected in the final state the result 
becomes 



s (i) 

9l,q 



4tt 



S K 



e/2 
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The timelike spin splitting function is given by 

A n , : , =c F [8( T i 7 ) + 

and the non-pole terms are given by 



- 4 - 4z + 65(1 - z) 



ln(l - z) 
1 - z 



i - z ; + 



l+z 2 

-2(1 + z) ln(l -z)+ 4 ln(z) + l-z 

1 — z 



(3.11) 



+<*(!-*)( -9 + 8C(2) 



(3.12) 



a (1) 

"1,9 



+2 



r/ln 2 (l - z) 



l + z 2 



§(^) + -^-«m>(i^ 



ln(z) ln(l-z) + ln 2 (z 



>(«) 



1 - z 

+(1 + z) (?C(2) - ^ ln 2 (l - z)) + (1 - z) (i ln(l 
+ ln(z))-^ + ^ + ( 5(l-z)(9-|c(2) 



1 - z 
-1 

^2 



- z 



(3.13) 



Notice that we also have to include in order a s terms proportional to e, which 
are represented by b^J Eq. (|3.1| ) and a^j. (k = q, g) since they contribute via 
mass factorization to the non-logarithmic parts of the coefficient functions. 
When the gluon is detected in reaction ( |3.9|) we get 



9i 



(i) 



at* 



'Q 



1 



m 



4tt \V , 

Here the timelike spin splitting function is given by 



8 -4z 



and the remaining coefficients are 



r (1) - 



[8 - 4z)(ln(l - z) + 21n(z)) - 16 + 12z 



(3.14) 



(3.15) 



(3.16) 
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l l,g 



C h 



(2-z)( ln 2 (l - z) + 4 111(2) ln(l - z) + 4 ln 2 (;z) 



+ (-8 + 6z)(ln(l - z) + 21n(z) - 2 J - (6 - 3*)C(2)] . (3.17) 

In order a 2 we have the following subprocesses. The non-singlet partonic 
fragmentation function receives contributions from the two-loop corrections 
to reaction (|3.7|) (fig. 5 of |§), the one- loop correction to reaction (|3.9| ) (fig. 
6 of H) and the double gluon bremsstrahlung process (fig. 7 of @) 



V -> q + q + g + g 



(3.18) 



where either the quark or the anti-quark is detected in the final state. Further 
we have to add the contribution from the partonic subprocess 



V-^q + q + q + q. 



(3.19) 



with the condition that the two quarks or the two anti-quarks which belong 
to the inclusive state are identical (see figs. 8,10 in ||). The result has the 
following form 



-NS,(2) 
9l,q 



An 
1 r 1 



'91 



r 1 



+^(A^«-A^«)-2/3b(^ + *g 



NS,(2) 
11 



+z 



The convolution symbol denoted by ® is defined by 



Cf® <?)(*) 



dzi / dz 2 5(z - z 1 z 2 )f(z 1 )g(z 2 ) ■ 



(3.20) 



(3.21) 



The terms c^l and a^ g already appeared in g[ l g (Eq. ( |3.10|) ). The order 



;(!) 



(1) 



non-singlet timelike spin splitting functions AP 1 ^ 1 -' and AP?^ 1 ' are 



11 



calculated in p8| (see also [E0]) and can be found in Eqs. (4.11) and (4.12) 
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of |J. Notice that they are equal to the spin averaged ones contrary to 
the singlet case which we will present below. The splitting function AP^q S '^ 
orginates from process (|3.19|) with two identical (anti-) quarks in the inclusive 
state. The same holds for the non-pole term c^' ( ' 2 ' ) ' ld []. Notice that c/i^'^ 
only appears if the amplitude G^ u (see Eq. [2.18|) ) is projected on the non- 
singlet part. If we project on the singlet part we get 



5 S,(2) 
9l,q 



~NS,(2) . ~PS,(2) 
9l,q + 9l,q 



(3.22) 



with 



-PS,(2) 
9l,q 



(i) 



+£AP#> ® ^1} + €' (2) + >P<?> ® afj + <?' (2) J , (3.23) 



where c lg and a lg are given in Eq. ( 13.141) . The purely singlet (PS) part 
originates from process ( |3.19| ) which proceeds by the exchange of two gluons 
in the t-channel (see figs. 8,9 in ||). Further z qq ^ = z^ q '^ + nfZ qq '^ where 
n f z qq originates from the light by light scattering graphs in fig. 1 in 



In the second order expression above the lowest order splitting function AP^ 
appears for the first time and it reads 



AP (o) = T f (-8 + 16z 
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(3.24) 



From our calculations and the expression given in Eq. ( |3.23 ) one can also 
infer the splitting function 

ap ps,(i) = CpTf f( 16 + 16z j ln 2 z _ ( 48 + n 2z ) \ nz _ 176 + I76z] . (3.25) 



The expression above agrees with the recent result obtained in ||20|| . The 
complete singlet spin timelike splitting function is then given by 



A pS,(l) _ A pNS,(l) 
qq qq 



AP 9 f « + AP^ 



qq 



(3.26) 



1 For the explicit expressions of 



NS,(2),nid , _NS,(2).id „ . orA A „n . 

and c l q see Eqs. (4. A. 32) - (4. A. 37) m 
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If the gluon is detected in the final state the one-loop corrections to reaction 
(|3.9|) and the contribution due to the gluon bremsstrahlung process in Eq. 
(|3.18p lead to the answer. 



-(2) 

9i,g 



'Air 



+2AP(J) ® (a« + 6W 



(3.27) 



In the second order expression above one encounters for the first time the 
lowest order splitting function 

a^ 0) = ca ( 8 (r^) + + 8 - 162 + y 5{1 " z) ) " l n f T ^ - z) - (3 - 28) 

From our calculation we also infer the second order spin timelike splitting 
function 



AP (1) 

99 



cl 



128 - 64z) (Li 2 (l - 2) + - In z ln(l - z) 



-— \n 2 z + -ln 2 (l -2) -C(2)l + (64-36z)lnz 
16 8 7 

-32(1 - 2) ln(l - 2) + 76 - 56z 

+C A C F [(32 + 16z)(Li 2 (-z) + lnzln(l + zj) + (-128 + 64z) 
x (Li 2 (l - z) + j In ^ ln(l -z) + - ln 2 (l - z)) - (48 + 8z) In 2 z 



(3.29) 



+ (128 - 48z)C(2) + (-32 + 88z) \nz + 32(1 - z) ln(l - z) 



+40 - 32,2 
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which agrees with the recent result in |20 . 

To obtain the finite coefficient functions in the MS-scheme one has to per- 
form mass factorization. The relations between the partonic fragmentation 
functions and the coefficient functions are then given by 



7 NS 
J qq 



-NS 
<?<? 



ACZ 



9u 



9i, 



2l®r M ®AC lj , + 7i / r (B ,<8)ACi, 



2Z, 



qq 



gq 



9 • 



9 ■ 



(3.30) 
(3.31) 
(3.32) 



Here the renormalization factor Z r is needed to restore the anticommutativ- 
ity of the 75-matrix on the level of the coefficient functions so that they are 
presented in the genuine MS-scheme. If we expand the transition functions 
in the unrenormalised coupling constant a s they read as follows 



-NS 

qq 



2W2 



47T 



e qq 



X 



I{iAP(>AP(;)-/3 AP(°)} 



2\ e 



+7{^( AP » S ' (1, - AP * ,(1 ')} 



6^2 

The singlet transition function is given by 



r: 



^NS ^PS 

T qq + 2n f T n 



qq - qq 1 ~' 4 / 1 qq > 

where the purely singlet (PS) part can be written as 



-PS 

qq 



^2 S 2 

■An 6 



M_ 

J 2 



2\ £ 



1 rl 



e 2 U 



{^AP(^Ap(o)} + i{iAP^«}]. (3.35) 



Finally we have 



99 



Air 



S, 





i'" 


-APW" 


U 2 , 







+( £ )8fl J 



(3.33) 



(3.34) 



18 



£ ^2 

We only need the following transition functions up to order a s 



(3.36) 



r -*±S f^- 
9ff ~47r £ U 2 



2\ e/2 



2£ 99 



(3.37) 



M 



2\ e/2 



Tgg — 1 + A S e ( „ 



47T 



99 



(3.38) 



Further we have to perform coupling constant renormalization which is also 
carried out in the MS-scheme. Up to order a 2 it is sufficient to replace the 
bare coupling constant by 



a s _ a s (M 2 ) s | a s (M 2 ) 2(3 0s ( M^\ e/ \ 
4ir 4n ^ An e 6 \ ji 2 j ' ' 



(3.39) 



The coefficient functions (MS-scheme) have the following representations 



M 47T 



■ + {Wf" - AP,f <") - A,c« 



+ iAP<;)®c«:»}L M + e^' 2 »'" M -^ 



_NS,(2),id 
1 



(3.40) 



PS 



-2 
J M 



+{iAP-' (1) + lAP4«)»cS}LM + ct S « 



(3.41) 
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4n 



7&y 

,9 



a 



s \2 



(i) 



^( 2 ) 
9 



with the definitions 



(3.42) 



(3.43) 



In the above expressions M denotes the factorization scale which has been put 
equal to the renormalization scale defined by R. In the case these scales are 
chosen to be different the corresponding coefficient functions can be derived 
from expressions fl3.41 )-( j3.43[ ) by the replacement 



a s {M 2 ) = a s (R 2 



1 , <*,(&) a . R 2 



As has been indicated in Eq. (|2.21j) AC^ q (|3.40| ) is equal to Cq " and the latter 



(3.44) 



is calculated in || (see also JK| ). Since z^f^ is not known yet, except for the 

first moment (see below), we can only present c^ S q '^ + z^i 5 ^ in Cf q (see Eqs. 

( |3.41|) , (|A.1| ) ) but not ^ S q and z^'^ separately. Unfortunately we cannot 
infer the latter from a ratio of two partonic fragmentation functions like we 
could do for z^ g s '^ in Eq. (|3.4| ). In the purely singlet case one has to compute 
the three-loop operator matrix elements for general moments (spin) given by 
the light by light scattering graphs in fig. 1 of [25]. The renormalization 



of these graphs is ruled by the Adler-Bardeen theorem [27] which is hard 
to impose on the level of the second order partonic fragmentation functions. 
For that one has to calculate the latter quantities in third order which is a 
tantalizing enterprise. Fortunately this does not hamper the computation of 
AC lg (|3.42 ) up to second order, in which z qq '^ does not show up, and the 
result can be found in Eq. ( |A.2| ). 

Finally we want to mention another problem which concerns the renor- 



malization constant Z 



NS, S) in Eq. ( |3.1| ) which is needed to restore 
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the Ward-identities broken by the HVBM-prescription for the 75-matrix. Af- 
ter coupling constant renormalization it reads 



qq 



5(1 



,e/2 



z$t(z)+eb<»(z) 



"Air' 



(3.45) 



where and byj are given in Eqs. 



qq 



and ( P-6| ) respectively. In the 



non-singlet case z^ q > can be inferred from Eq. ( |3.4j ) and it is presented in 



Eq. ( |A.3j ) . The first moment of Eq. (|3.45 ) is equal to 



7 NS 
J qq 



a: 



1 + -^S e C F 

4:71 



4+ 13e 



+ 



, a s 2 
"4^ 



Cf.{12 + 16C(2)} 



441 

~3~£ 



107 161 4 

\ + nfCpTfi — \ 

9 1 t tX 3 £ 9 J 



(3.46) 



where we have put M — \i. Comparing the result above with the one obtained 
in 



1 1 ] given by 



7NS 

"99 



1 + £-s g c F 

47T 



-4 + 5e 



CM 22} 



_ _ . 441 107, „ ^ .161 4, 
H-C,C f {- T ---} + .vC F T /{T - + 5 } 



(3.47) 



we observe a discrepancy which shows up in the order a 2 s coefficient of the 
colour factor C\. Notice that the result in Eqs. (8), (11) of [11 is ob- 
tained by computing the ratio between the vector and the axial-vector form 
factor where it was needed to restore the Ward-identity for the non-singlet 
axial vector current broken by the HVBM-prescription. A discrepancy is 
also oberved when Eq. (|3.4j) is computed for deep inelastic scattering where 



q = —Q (see ||). Here one gets an answer which differs from the result in 



Eq. (|3.46p but agrees with Eq. ( J3.47[ ) . From the above we conclude that the 



renormalization constant Z^„ is process dependent except for leading order. 



qq 



The difference of this constant between deep inelastic structure functions 
(spacelike process) and the fragmentation functions (timelike process) is of 
the same origin as the one discovered for spacelike and timelike splitting func- 
tions in ||28|| . The process dependence of the renormalization constant above 
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is in contrast to what is known about the usual renormalization constants 
computed in cases where the regularization procedure does not violate the 
Ward-identities. In |25[ also the first moment of Z qq has been computed. 



The difference between Z^ q and Z qq only shows up in the nfCpTf-paxt of 



Eq. ( |3.47D . In the singlet case one has to add to the latter expression the 
first moment of rifZ qq '( 2 > which equals 3nfCpTf. Since the leading terms in 
Z qq (r = NS, S) are process independent we can assume that this also holds 

for z qq '^ so that at least we know the first moment for AC^ q (Eq. (|3.41|) ). 

Summarizing the above we have computed the order a 2 contributions to 
the coefficient functions corresponding to the spin fragmentation functions 
g^(x, Q 2 ) (i = 1,4,5) measured in electron-positron annihilation in which 
the hadron H is polarized. In addition we obtained the second order split- 
ting functions AP P ? S, W, AP^ which agree with the recent results published 
in . Unfortunately we could not get the remaining ones given by and 
which could be calculated by the cut vertex method in PC|| . This is be- 
cause the splitting functions P^ pw always appear one order higher in the 
coefficient functions than P qq , Pjjq- This property can be traced back to the 
fact that the electroweak vector bosons never directly couple to the gluon. 
Finally we want to mention that for a complete next-to-next-to-leading order 
analysis one also needs the three-loop splitting functions which have to be 
combined with the coefficient functions in Eqs. ( |3.40|) - (|3.42|) . The computa- 
tion of these splitting functions for spacelike and timelike processes is one of 
the outstanding enterprises which have to be still carried out. 
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Appendix A 



In this section we only present the coefficient functions AC X and AC 1)9 since 
the non-singlet ones, given by AC 



77NS 



= 1,4,5), are already known in the 
literature via the relations in Eq. ( |2.21| ). The purely singlet coefficient equals 



n f {— fC F T f [{8(1 + z) lnz + 20(1 - z))l 
An - 1 > 



2 

M 



+ {l6(l + z) (li 2 (l - 2) + In 2 ln(l - 2)) + 24(1 + z) In 
+ (8 - 88z) In ^ + 40(1 - z) ln(l - z) - 136(1 - z)}l a/ 
+16(1 + z)(3Si )2 (l - z) - Li 3 (l - z) + 31nzLi 2 (l - z) 



1 3 
+ ln(l - z)Li 2 (l - z) + - In z ln 2 (l - z) + - In 2 z ln(l - z) 

+ ^ In 3 z - C(2) lnz) + (8 - 88z) (Li 2 (l - z) + In z ln(l - z 



32 32 
- (— + 32 + 32z + — z 2 ) (Li 2 (-z) + In z ln(l + z)) + ( - 26 



16 



58z + — z 2 ) In 2 z + 20(1 - z) ln 2 (l - z) + ( - 72 + 40z 



with 



32 



3 
472 



- 2 C(2) 



364 356 



z lnz- 136(1 - z) ln(l - z) 



(1---) -n f Qrz^( Z) 



(A.l) 



^( 2 )( z ) = C F T / [-4(l + 2z)ln 2 z + 8(-l + 3z)lnz + 16(l-z)l . (A.2) 



The gluonic coefficient function becomes 



a 



s \2 



"A.-K' 



Cf{ ((-8 + ln z + (16 - 8z) ln(l 



- z 
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+Qz)L 2 M + ((48 - 24z)Li 2 (l - z) + (-24 + I2z) In 2 z + (48 
-2Az) \nz\n(l-z) + (32 - 16z) ln 2 (l -2) + (-64 + 32z)C(2) 
+ (128 + 4z) lnz + (-80 + 682) ln(l - z) + 132 - 148z)l m 

/ 96 \ /64 \ 

+ ( + 64 - 96z)U 3 (-z) + (— + 128 + 64z)Si, 2 (-z) 

Z ' Z 

+ (— - 240 + 120z)Si ;2 (l -*) + (— + 128 + 64z) 

Z Z 

x ln(l + z)U 2 (-z) + (-48 + 2Az) lnzLi 2 (l - z) + (32 - 16z) 

32 

x ln(l - z)U 2 (\ - z) + ( 64 + 32z) \nzU 2 (-z) + (8 - 4z) 

16 

x In 2 z ln(l - z) - (— + 32 + 16z) In 2 z ln(l + z) + (40 - 20z) 

x lnzln 2 (l -z) + (— + 64 + 32z) lnzln 2 (l + z) + (- ^ 

v z / v 3 

+f z) In 3 * + (f - f *) ln 3 (l - z) + (16 - 8z)C(2) lnz 
+ ( - — + 48 - 24z)((2) ln(l - 2) + (— + 64 + 32z)C(2) 

Z ' z 

32 

x ln(l + z) + ( + 80 - 104zV(3) + (-32 + 20z)Li 2 (l - z) 

^ z ' 

+ (-64 + 52z) In z ln(l - 2) - f — + 64 + ^z 2 ) (Li 2 (-z) 

v 3z 3 / v 

32 

+ In z ln(l + z)) + (128 - 43* + y z 2 ) In 2 z + (-72 + 54z) 
x ln 2 (l - z) + (32 - 60* - y, 2 )C(2) + (y + ^) In* 

24 



+ (^ + 88 - 108z) ln(l - *) + - 
^ z ' 3 3 ' 

+C A C F { ( - (32 + 8z) In z + (16 - 8z) ln(l - z) - 48(1 - z)) 
x4 + ( - 96Li 2 (l - z) - 32(1 + z) lnzln(l - z) + (32 + 16z) 
x(Li 2 (-z) + lnzln(l + z)) - (112 + 24z) In 2 z + (16 - 8z) 

x ln 2 (l - z) + (96 - 32z)C(2) + (-32 + 136z) lnz + (-128 

1 -z N 



+112z) ln(l - z) + 232 - 224z)l m + (32 + 16z)(Li 3 (--^—^j 
-Li 3 (^|) ) + (80 + 40z)Li 3 (l - z) - + 320) 



xSi )2 (l - z) + (— + 32 + 16z)Li 3 (-z) - (— + 64 + 32z) 

Z Z 

xS li2 (-z) - (288 + 16z) lnzLi 2 (l - z) - (80 + 8z) ln(l - z) 

xLi 2 (l - z) + (32 + 16z) ln(l - z)Li 2 (-z) + ( + 32 + 16z) 

^ z ' 

32 16 
x In zLi 2 (-z) - (— + 64 + 32z) ln(l + z)Li 2 (-z) - (— + 32 

Z ' z 

248 

+ 16z) lnzln 2 (l + z) - (96 + 32z) ln 2 zln(l-z) - (— 

44 \ 

+— z) In 3 z + (32 + 16z) In z ln(l - z) ln(l + z) - (24 + 12z) 

x In z ln 2 (l - z) + (- + 48 + 24z) In 2 z ln(l + z) + (- - -z) 

^z ' ^3 3' 

16 



x 



ln 3 (l - z) + (256 - 48z)C(2) lnz + (— + 16z)C(2) ln(l - z) 
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— + 32 + 16z)C(2) ln(l + z) + (— + 24 + 20zV(3) 



1 fi 

+ (64 + 56z)Li 2 (l - z) + (72 + 38z - — z 2 ) In 2 z + 



176 32 
(-56 + 48z) ln 2 (l - z) + (— + 64 + 32z + — z 2 ) (Li 2 (-z 

+ In z ln(l + z)) + (-48 + 136z) In z ln(l - z) + (48 - 16z 



+^ 2 )C(2) + ( 



844 128 



3 3 

, m \ 124 76 1 
x ln(l - z )- — + -j z \ 



:) In : + ! - - + 276 - 240z 



(A.3) 



The renormalization constant in Eq. (|3.46|) needed to restore the Ward- 
identities, which are broken by the HVBM prescription of the 75-matrix, is 
equal to 



Z N „ s (z) = 5(1 - z ) + r±C F -8(l-z) +e{-8(l-z)lnz 



a 



99 



47T 



-4(1 -z) ln(l-z) + 10-8z)}" 

+ (^) 2 [C 2 F { - 16(1 - z) - (8 + 16z) lnz + 16(1 - z) In 2 z 
-16(1 -z) lnzln(l-z) 

+C^C F { - (1 - z) - ^(1 - z) + 8(1 - z) C(2) 



+ (-y + ~z)lnz -4(1 - z) ln 2 zj 
+n f C F T f {- £ f(l - z) + I (1 - z) lnz + ^(1 - z)} 

" ^CUC F ){8(1 + z)(4Li 2 (-z) + 41nzln(l + z) + 2((2) 



26 



In 2 



3 In z 



56(1 -z)] 



(A.4) 



Where the last part originates from the process in Eq. ( |3.19| ) with two 
identical particles in the inclusive final state. In the expressions above the 
definitions of the Riemann zeta-functions ((n) and the polylogarithms Li n (z), 
Sn,m{z) can be found in 
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